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ABSTRACT 

Radio emission, polarization, and Faraday rotation maps of the radio jet of the 
galaxy 3C 303 have shown that one knot of this jet carries a galactic-scale electric 
current and that it is magnetically dominated. We develop the theory of magnetically 
dominated or Poynting-flux jets by making an analogy of a Poynting jet with a trans- 
mission line or waveguide carrying a net current and having a potential drop across 
it (from the jet's axis to its radius) and a definite impedance which we derive. The 
electromagnetic energy flow in the jet is the jet impedance times the square of the jet 
current. The observed current in 3C 303 can be used to calculate the electromagnetic 
energy flow in this magnetically dominated jet. Time-dependent but not necessarily 
small perturbations of a Poynting-flux jet are described by the "telegrapher's equa- 
tions." These predict the propagation speed of disturbances and the effective wave 
impedance for forward and backward propagating wave components. A localized dis- 
turbance of a Poynting jet gives rise to localized dissipation in the jet which may 
explain the enhanced synchrotron radiation in the knots of the 3C 303 jet, and also 
in the apparently stationary knot HST-1 in the jet near the nucleus of the nearby 
galaxy M87. For a relativistic Poynting jet on parsec scales, the reflected voltage wave 
from an inductive termination or load can lead to a backward propagating wave which 
breaks down the magnetic insulation of the jet giving |E|/|B| > 1. At the threshold 
for breakdown, |E|/|B| = 1, positive and negative particles are directly accelerated in 
the E X B direction which is approximately along the jet axis. Acceleration can occur 
up to Lorentz factors ~ 10^. This particle acceleration mechanism is distinct from 
that in shock waves and that in magnetic field reconnection. 

Key words: galaxies: jets — accretion discs — magnetic fields — acceleration of 
particles 



1 INTRODUCTION 

A fundamental open question of astrophysical jet models is 
how energy extracted from the accretion flow close to the 
black hole event horizon is transported along narrow jets to 
much larger distances. The total energy carried by the jets of 
active galaxies is estimated to be a non- negligible fraction of 
the massive black hole formation energy, ~ O.lMhhC^ (Kron- 
berg et al. 2001). The jets are initially highly relativistic and 
low-density and for this reason are magnetically dominated 
or force-free. In this limit a negligible fraction of the power 
is carried by the particle kinetic energy. That is, the energy 
outflow from the accretion disc is in the form of a coUimated 
'Poynting-flux jet' as proposed by Lovelace (1976) and sub- 
sequently studied in many papers (Benford 1978; Lovelace, 



Wang, & Sulkanen 1987; Lynden-Bell 1996; Li, et al. 2001; 
Lovelace, et al. 2002; Lovelace & Romanova 2003; Nakamura 
et al. 2008). 

An illustrative model for the formation of a relativistic 
Poynting jet is sketched in Fig. 1. A large scale magnetic 
field is assumed to thread the accretion disc. As a result 
of the differential rotation of the disc the initial field loops 
opens up giving a current outflow (or inflow) / along the 
spine of the jet of cylindrical radius rj initially of the or- 
der of the Schwarzschild radius of the black hole (Lovelace & 
Romanova 2003). A three-dimensional view of a numerically 
calculated magnetic field is shown in Fig. 6 of Lovelace et al. 
(2002). In this model the jet power comes from the accretion 
disc. The associated toroidal magnetic field is responsi- 
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ble for coUimating the jet. An equal but opposite "return 
current" flows inward (or outward) at much larger radial dis- 
tances from the jet axis so that the net current outflow from 
the source is zero. Because the jet current and the return cur- 
rent have opposite signs they repel as a result of their mag- 
netic interaction mediated by the toroidal magnetic fleld. 
This repulsion between the jet and its return current has 
been demonstrated in magnetohydrodynamic (MHD) simu- 
lations (Ustyugova et al. 2000; Nakamura et al. 2008). The 
energy flow in the jet is unidirectional and is carried predom- 
inantly by the Poynting flux z ■ (E x B)//io — Er x B^/^o 
(MKS units), where Er is the radial electric field of the jet. 
The power in the Poynting jet may come predominantly 
from the spin-down of the black hole (Blandford & Zna- 
jek 1977). This process, which depends on the presence of 
a magnetized accretion disc around the rotating black hole, 
has been studied with general relativisitic MHD simulations 
(e.g., McKinney 2006; Beckwith, Hawley, & Krolik 2008; Tc- 
chekhovskoy, Narayan, & McKinney 2011). 

Rotation measure gradients observed on parsec scales 
close to the nuclear central black hole have given evidence 
of electric current flow along the jets (Asada et al. 2002; 
Gabuzda, Murray, & Cronin 2004; Zavala & Taylor 2005). 
More recently, on a scale 10* times larger, radio emission, 
polarization. X-ray, and Faraday rotation maps of the radio 
jet of the galaxy 3C 303 have shown that one knot of this 
jet has a galactic-scale current of ~ 3 x 10^* Ampere flowing 
along the jet axis (Kronberg et al. 2011). The physical pa- 
rameters of this knot derived from the observations indicate 
that the jet is magnetically dominated, that is, a Poynting 
jet. The existence of an axial current in this very large scale 
plasma flow suggests an electrical circuit analogy for the jet. 

A further possible indication of a Poynting jet is the 
jet near the nucleus of the galaxy M87 which has been ob- 
served to have an apparently stationary but strongly time- 
dependent knot HST-1 in its jet (Biretta et al. 1999). This 
knot showed a dramatic peak in emission in about 2005.2 in 
the radio, optical, and X-ray bands (Harris et al. 2009). 

In Sec. 2 we model a Poynting jet as a transmission 
line carrying a net axial current and having a potential 
drop across it. Further, we derive the transmission line 
impedance and the electromagnetic energy flow in the jet. 
Time-dependent but not necessarily small perturbations of 
a Poynting-flux jet (§2.2) are described by the "telegraphers' 
equations" (Heaviside 1893). These are wave equations for 
the current and voltage across the line. The voltage and cur- 
rent consist in general of forward and backward propagating 
components. A localized irregularity of a Poynting jet (§2.3) 
can give rise to localized dissipation in the jet which may 
explain the enhanced synchrotron radiation in the knots of 
the 3C 303 jet and in the apparently stationary knot HST-1 
in M87. 

Sees. 2.4 & 2.5 consider highly relativistic Poynting jets 
on parsec-scales jets such as HST-1. In this case the reflected 
voltage wave from an inductive load can lead to a backward 
propagating wave which breaks down the magnetic insula- 
tion of the jet, giving |E|/|B| ^ 1. The threshold for this 
breakdown is |E|/|B| = 1. At this threshold, positive and 
negative particles are directly accelerated in the E x B di- 
rection which is approximately along the jet axis and can be 
accelerated to very high Lorentz factors. Sec. 3 gives conclu- 
sions of this work. 




Figure 1. Sketch of the magnetic field configuration of a Poynt- 
ing jet adapted from Lovelace and Romanova (2003). The bot- 
tom part of the figure shows the initial dipole-like magnetic field 
threading the disc which rotates at the angular rate Q{r). 

2 THEORY 

In cylindrical (r, (j), z) coordinates with axisymmetry as- 
sumed, the magnetic field has the form B = Bp -I- B^<j> , 
with Bp = BrT + B^z, and Br = ~{l/r){d'i' /dz), and 
B^ = {l/r){d^ /dr). Here, *(r, z) = rA^{r,z) is the 
flux function. A simple form of this function is ^'(r, 0) = 
(l/2)r^_Bo/[l -|- 2(r/ro)^], where -Bo is the axial magnetic 
fleld strength in the center of the disc, and ro is the ra- 
dius of the O— point of the magnetic fleld in the plane of 
the disc as indicated in Fig. 1. This is taken to apply 
for r > even though it is not valid near the horizon of 
the black hole. The contribution from this region is neglig- 
ble for the considered conditions where (rg/ro)^ <C 1, with 
Tg = GM/c? . For a corotating disc around a Kerr black hole 
the disc's angular velocity viewed from a large distance is 
n = [c^ /{GM)]/[a, + (r/rg)^/^], for r > r^s where r^s is 
the innermost stable circular orbit and a, is the spin pa- 
rameter of the black hole with ^ a, < 1. 

At large distances from the disc (z 3> ro) the flux func- 
tion solution of the 'force-free' Grad-Shafranov equation is 



[27^ (r2 - 1) ]2/3 ' 

(Lovelace & Romanova 2003). The 'force-free' solution is 
applicable for conditions where the kinetic energy-density of 
the plasma is much less than the magnetic energy-density. 
Here F is the Lorentz factor of the jet, f = r/ro, ^ = ^7^*0 
with 'I'o = ro-Bo/2, and TZ = ro/vg. This dependence holds 
for fi = [2(r2 - l)]i/V7^ < r < fa = [27^(^2 - l)]i/V3^'''*- 
At the inner radius fi, — 1/TZ^, which corresponds to 
the streamline which passes through the disc at a distance 
r = Tg. For f < fi, we assume ^ cc , which corresponds 
to Bz = const. At the outer radius fa, ^ — (*I')max ~ 1/3 
which corresponds to the streamline which goes through the 
disc near the O— point at r = ro. Note that there is an 
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appreciable range of radii if TZ^^^ ^ 1. Also note that for 
z ^ ro, radiation drag on the lepton component of the jet 
is negligible. 

For fi < r < f2, the field components of the Poynting 
jet are 



Er = -\/2 (r^ - 1) 



1/2 



= -V2 r B^ 



and 



p-2/3 



3 [27^(^2 - 1)]2/3 • 

This electromagnetic field satisfies the radial force bal- 
ance equation, dB^/dr + {l/r^)d[r^{Bl - Er)]/dr = 0, 
for an a^xisymmetric, translationally-symmetric, and time- 
independent force-free field. Figure 2 shows the field profiles. 

Appendix A describes a different force-free jet model 
with electric and magnetic fields similar to those in a com- 
mon coaxial transmission line. 

At the jet radius r2, there is a boundary layer where 
the axial magnetic field changes from B^ (r2 — e) to zero at 
r2 + e, where e ^ r2 is the half-width of this layer. The 
electric field changes from Er{r2 — e) to zero at r2 -I- e. 
The toroidal magnetic field changes from B^{r2 — e) to 
B^{r2 + e) where this change is fixed by the radial force 
balance. Thus for r > r2, we have Er = Q, Bz = Q, and 
B4, = VSB,{r2 - e)(r2/r) = ^/(3/2)r-'' B^^^ri - s){r2/r). 
Equivalently, B^{r2 + s)/B^{r2 - s) = ^/iJ^jV. 

The toroidal magnetic field for r > applies out to 
an 'outer radius' rs where the magnetic pressure of the 
jet's toroidal magnetic field, B'^{r[i)/STv = pcx, balances 
the external ram pressure Pcx ~ Pcx + Pc^iidr-.i/ dt)^ , where 
Pex = riex/cBrex is the kinetic pressure of the external inter- 
galactic plasma and pex is its density. The outward propa- 
gation of the jet will be accompanied by the non-relativistic 
expansion of the outer radius, drs/dt > 0. 

We take as the 'jet current' the axial current /o flow- 
ing along the jet core r < r2 — e. From Ampere's law, 
B^{r2) = — 27o/(c r2) or in convenient units, S<^[G] = 
— /o[A]/(5r[cm]). The net current carried by the jet {r ^ 
r-2 +e) is /not = ys/Ir-iJo. 

Using equations (2), the energy flux carried by the 
Poynting jet can be expressed as 



Ej 



Jo 



rdrErB^ = ZqIq 



c 2n 
2' Ho. 

where the square bracket is for cgs or MKS units, and 



Zo = [cgs] = 90 M - ^ 



1/2 



n [MKS] 



(3) 



(4) 



where = U,/c = (1 - r-^)^/^ Here, Zq is the DC 
impedance of the Poynting jet. The conversion to MKS units 
is — > (4:TT)~^(po/€oy^^ = son. Earlier, the impedance 
of a relativistic Poynting jet was estimated to be ~ 
(Lovelace 1976). 

We can apply these concepts to the jet in 3C 303 where 
the observed axial current in the E3 knot is 3.3 x 10^^ A 
(Kronberg et al. 2011), Thus the electromagnetic energy flux 
is Ej Ri 8 X 10*^/3 erg s~^. This energy flux is much larger 
than the photon luminosity of the jet of 3.7 x 10"*^ erg s~^ 
integrated over 10® to lO^'' Hz (Kronberg et al. 2011) as- 
suming /3 = Uz/c is not much smaller than unity. For the 
E3 knot the jet radius is r2 w 0.5 kpc so that B^(r2) w 0.4 




Figure 2. Radial dependence of the field components of a Poynt- 
ing jet from equation (8). The inner jet radius is assumed to be 
ri = 0.05r2, the boundary layer thickness is £ = 0.01r2, and 
r = 5. 



mG. The E3 knot is about 17.7 kpc in projection from the 
galaxy nucleus. 

2.1 Transmission Line Analogy 

Here, we interpret the Poynting flux jet in terms of a trans- 
mission line as proposed earlier by Lovelace & Ruchti (1983). 
The different physical quantities are measured in the 'labo- 
ratory' frame which is the rest frame of the plasma outside 
of the jet at r ^ r2. The effective potential drop across the 
transmission line is taken to be 



1 

Vo = --ro J dfEr{f) 



(5) 



ro Bp 
3V4 ' 

where the factor of one-half accounts for the fact that the 
transmission line does not consist of two conduction sur- 
faces. 

The axial current flow of the jet is 



T 1 R ^ ^ ^0 



(6) 



with Zo given by equation (4). The units of equations (4) and 

(5) are cgs. In MKS units note that a current Jo = 3 x 10^® 
A gives a voltage Vo = 2.7 x 10'^"/? V. The energy of ions 
accelerated across this voltage is large enough to account 
for ultra high energy cosmic rays (Lovelace 1976; Biermann, 
Kang, & Ryu 2001; Ostrowski 2002). 

There is a corresponding electric field energy per unit 
length of the jet which in MKS units is 

47reo 



We 



\drEl = ^CVo 



Here, C = 



(7) 



is the capacitance per unit length in Farads per meter and 
eo = 8.854 x 10"^^ F/m. 

The magnetic energy per unit length of the jet in MKS 
units is 

= JL rrdriBl + B^z) + - rrdrBlJny , 

fJ-O Jo tM> Jrn \ r J 
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(8) 



Here, B</>+ is the toroidal field at r2 + e. Carrying out the 
integrals we find 



L = 



3^0 

4-K 



(9) 



which is the inductance per unit length in Henrys per meter 
with /io = 47r X 10~^ H/m. 



z 






i 






Vo 


AV_-<-| 


z^- I^^AV, 



(10) 



2.2 Telegraphers' Equations 

Time and space (z— )dependent perturbations (not necessar- 
ily small) of a Poynting-fiux jet are described by the Telg 
rapher's equations, 

dAV _ 1 QAJ dAI _ 1 aAV 

dt ^ ~c dz ' ~ar ^ ~I 

where (AV, A/) represent deviations from the equilibrium 
values {Vo, lo) (Heaviside 1893; Bergeron 1977; Samokhin 
2010). The equations can be combined to give the wave equa- 
tions 

where 



{AV, A/) = 



(11) 



1 



1 + 



1 

2r2 



1 

2r2 



In 



-1/2 



(12) 



is the phase velocity of the perturbation. The general solu- 
tion of equation (10) is 



AV 
A7 

with 

Here, 

Z = 



AV+{z ~ u^t) + AV-{z + u^t) , 
AI+{z - u^t) + AJ_(2 + u^t) , 



(13) 



AV+ 



= 90 



2AI+ and AV_ 



-ZAI- 



1 

2f2 



1 

2r2 



1/2 



J7 [MKS], (14) 



is the wave impedance of the jet. 



2.3 Irregularities in the Transmission Line: 

A localized irregularity may appear in the transmission line 
due possibly to an instability at t > as diagrammed in Fig. 
3. The irregularity can be modeled as an extra impedance 
Zi or 'load' across the transmission line at z = 0. This 
impedance is considered to go from Zt{t < 0) = oo to a, 
constant value Ze for t > 0. In general Zt is complex with 
an imaginary component (reactance) which can be positive 
(capacitive) or negative (inductive). On either side of the 
discontinuity the line is assumed to have the impedance Z 
given by equation (14). On the upstream side of Ze (^ < 0), 
the line voltage is Vo + AV- , where AV- is the backward 
propagating wave. The current on this part of the line is 
7o -I- A/_. There is no forward propagating wave for the 
conditions considered here. On the downstream side of Ze, 
the line voltage is Vo + AVt and the current is 7o + A7t, where 
(AVt, A7t = AVt/Z) represent the transmitted wave. 

The standard conditions on the potential and current 



Figure 3. Diagram of an irregularity in the transmission line 
represented by the impedance Ze{t) where Zi{t < 0) = oo and 
Ze{t > 0) is finite. Here, V_ the wave reflected off of Ze, Vt the 
transmitted wave, and Vq is the line voltage for t < 0. 



Figure 4. Diagram of a transmission line of impedance Z termi- 
nated by a load impedance Zj^ with an inductive component. 



flow at Ze {z = 0) give Vo + AV- = Vo-f- AVt and 7o + A7_ = 
Ie + Io + Alt , where le — ( Vo -|- A V- )/.E< is the current flow 
through Ze- In this way we find 

AV_ = ^^Vo = AV. . (15) 

Note that for Ze — >■ oo, both AV_ and AVt tend to zero. 
The power loss rate in the load Ze is 



{Vo + AVtf 4Ze 



Ze 



(16) 



(22, +2)2 " • 

We assume that this power goes into accelerating charged 
particles which in turn produce the observed synchrotron ra- 
diation. This power could account for the emission of the E3 
knot of 3C 303 (Kronberg et al. 2011) and, on parsec scales, 
the emission of the HST-1 knot in the M87 jet (Biretta et 
al. 1999). 

The power Ee is equal to the difference between the 
Poynting flux for z < 0, 



[Vo^ - {AV-f\ 



Z 



(17) 



S- = (V) + AV-)(7o-hA7_) = 
and the Poynting flux for 2: > 0, 

^ (Vo + AV)(7o + A70 ^ [V-o- + 2VoAV + (AV)-] ^(^3) 



Z 



Using equation (17) we flnd that S 



2.4 Magnetic Insulation 

The Poynting jet is analogous to a magnetically insulated 
transmission line. In laboratory experiments magnetically 
insulated transmission lines involving crossed electric and 
magnetic fields are commonly used to transport pulses of 
high energy content onto small targets (e.g., Shope et al. 
1978; Samokhin 2010). The condition for "magnetic insu- 
lation" of a planar gap of width d between a conducting 
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cathode and anode with potential drop Vb and initially filled 
with a uniform magnetic field Bo is Bq > (214 Vb + Vo)d^^, 
where K = m^c^je « 5.14 x lO'^ V (Lovelace & Ott 1974; 
Ron, MondelU, & Rostokcr 1973). 

In the astrophysical limit where Vb ^ Ve, the condition 
for magnetic insulation is |E| < |B| for the case where field 
components have the same radial dependence. If the power 
flow in the jet is unidirectiorial, as assumed in equations (2), 
then the condition for magnetic insulation is 



M 

IBI 



1/2 



< 1 or 



Vb 



l+2f^) 



1/2 



(19) 



,r2 + i/2y ZrIo 

where Zr = 90Q is a reference impedance. This inequality 
is always satisfied because the fields are time-independent 
and force- free (peE -|-J x B/c = 0). 

Consider now the case where the power flow is not 
unidirectional and where the fields are time and space- 
dependent, i.e., V{z,t) = AV++AV-, I{z,t) = AI++AI-. 
We assume > 1 so that {B^/B^f = ©(F'^) < 1. To 
a first approximation the electric and magnetic fields corre- 
spond to those in a common transverse electric and magnetic 
(TEM) mode in a coaxial cable. Unlike a cocixial cable where 
the two conducting surfaces are separated by a dielectric, 
the Poynting jet has distributed charge and current densi- 
ties. Because \Er/B^\ = V/{2I), the insulation condition 
lEI < IBI can then be written as 



(20) 



where C = 1 + 0{r ^)isa critical dimensionless number 
slightly larger than unity. 



2.5 Consequences of Breakdown of Magnetic 
Insulation 

Here we discuss conditions where the magnetic insulation 
breaks down in a localized region near the termination of a 
Poynting jet due to a load 2l which has an inductive com- 
ponent (Fig. 4.). An essential condition for this to occur is 
that the wave propagation not be unidirectional. The load 
is assumed to be stationary or sub-relativistic in the labo- 
ratory frame owing to the ram pressure from the external 
plasma. The breakdown of magnetic insulation in laboratory 
transmission lines is discussed by Gordeev (1978). 

The consequences of the breakdown can be dramatic: 
Nuclei can be accelerated to energies of the order of eV ~ 3x 
10'^" eV while leptons may be accelerated to smaller energies 
owing to radiative losses and interactions with radiation. 
This value of the potential comes from the jet parameters 
deduced for 3C 303 with Z = 90fi (Kronberg et al. 2011). 
(Note that magnetic insulation breakdown can also occur 
along the jet due to the refiected wave from a localized region 
of increased series inductance along the transmission line.) 

The voltage across the Poynting jet, V = AV+ + AV-, 
is assumed to be made up of an incident component AV+ 
propagatiiig along the jet channel with velocity «^ and a re- 
fiected component AV- with velocity — m^. The correspond- 
ing jet current is / = AI+ + A/_. The spatial profile of 
the incident wave AV+ is considered to be a smooth rise to 
a constant value (taken to be unity) as shown in the top 
panel of Fig. 5. Such jet onsets are observed in the form of 
the sporadic formation of parsec-scale components of radio 



V ^..^ 


AV. 




\ 



1 



- 



-10 


z 





V 

- zif\ 








ZM. 


-10 


z 






Figure 5. Reflection of an incident wave (AV^ , A/-|-) off of a sta- 
tionary load of impedance Zi^ and tfic reflected wave (AV_, A7_) 
for tfic circuit in Fig. 4. Tfie net voltage is V = AV+ -I- AV_ 
and tfie net current / = A7_)_ -|- A/_ . The top panel sfiows the 
voltages and the bottom panel the currents as well as the ratio 
V{ZI)~^. If this ratio is larger than a critical value, C (equa- 
tion 20), then the magnetic insulations will breakdown. For this 
example, 2^ = Z — ucC, where C is the inductance of the load. 



galaxies and quasars (e.g., Zensus et al. 1998). The onset 
may be explained by the global magnetic instability of the 
disc (Lovelace et al. 1994). 

In such cases the wave amplitudes can be written as 



AV±{z,t) 



-f 

J — c 



with AV±{ij:y 



dui AV±{to) exp[iuj{±z/u^ - t)] , (21) 

: AV±{—lli) and with 

AF+(w) (22) 



^^-^'^^ Z.^Z 

where .2^(0;) is in general frequency-dependent with real and 
imaginary parts. Figure 5 shows a case where Zl = R — iwC 
is due to a resistance R assumed equal to Z in series with 
and inductance C An estimate of C is simply r2L, where r2 
is the jet radius and L is given by equation (9). 

The top panel of Fig. 5 shows that there is a backward 
propagating positive voltage wave AV^. The rise-time of 
the wave is equal to that of the incident wave while the 
decay-time is equal to the L over R time of the load. The 
voltage wave is positive because the inductance of the load 
appears initially to the incident wave as a high impedance. 
At the same time the refiected current A/_ is negative. The 
net voltage and current arc such that there is a backward 
propagating wave with V/{ZI) > C. This corresponds to a 
breakdown of the magnetic insulation. 

The condition V/(ZI) > C can exist only transiently 
(for a time < r2 / c) because for Er < the insulation break- 
down will lead to a radial inflow of positive charge with speed 
c(l — B^/E^y^^ and an opposite radial outflow of negative 
charge. At the same time both positive and negative charges 
drift in the E x B direction with speed c\B/E\. This wiU 
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bring the transmission line to the threshold condition where 
V/{ZI) =C or |E| = |B|. 

Both positive and negative charged particles in the re- 
gion where |E| = |B| are accelerated in the E x B direction 
which is approximately the +z direction. The acceleration 
is highly efficient in that a particle exposed to the wave for 
a time At is accelerated to a Lorentz factor 

1 f&eBMV'^ 



7exb ~ -'- + 2 V 



(23) 



which is much larger than unity. This formula is easily de- 
rived from the equations of motion. The electromagnetic en- 
ergy in the backward propagating wave is dissipated into 
kinetic energy of the accelerated particles. 

As a numerical example relevant to the HST-1 compo- 
nent in M87, we consider r2 = 5 pc and B = 0.5 mG (Har- 
ris et al. 2009) and assume that the jet consists of electrons 
and positrons. De Young (2006) discusses different models of 
extra-galactic jets involving plasmas dominated by electrons 
and positrons, electrons and ions, and electromagnetic fields 
(i.e., Poynting fiux jets). For the electron-positron case, 

7EXB « 1.2 X 10^ i ( ^-^Y' , (24) 



5 X 10»s 



where the time At is normalized to r2/c « 5 x 10*s. The ac- 
celerated particles will in a time ~ r2/c run into the region 
of the load where the magnetic field has a spread of direc- 
tions. The L/R time ~ 5 x lO^^s is much longer than the 
dissipation time for the backward propagating wave. The 
synchrotron radiation peak from these high-7 leptons is at 
I'syn = SeB'f^ /2meC w 3 X lO^'^Hz or a photon energy w 1300 
eV. 

If the electromagnetic radiation of the source frad 
is entirely due to the accelerated particles, then frad ~ 
NiirieC^'YExB, where Ni is the rmmber of accelerated lep- 
tons. Estimating Ni ~ Trricn^ and taking frad = 10'*^erg/s 
(Harris et al. 2009) gives the lepton number density ne ~ 
4.5 X 10~® cm~^. Prom this we obtain the so called 'plasma 
beta', I3p = AnnfyExBrneC^ / B'^ ~ 2.2. This implies that the 
load region of the jet is not force-free and that it will ex- 
pand. Prior to the acceleration of the particles when the 
power flow in the jet is unidirectional, /3p = 1.9 x 10~^r so 
that the jet is dominantly force-free. Synchrotron and other 
losses will give rise to a power law distribution of lepton 
energies oc 7"' with q ^ 2 and with a cutoff at jExBrrieC^ . 

The estimated Lorentz factor (24) is larger than a crit- 
ical value denoted 7c set by the condition that the inverse 
Compton scattering of a lepton off a synchrotron photon 
(energy Esyn), which gives an inverse Compton photon (en- 
ergy eic w 7^esyn), be such that the subsequent scattering 
of the synchroton and secondary inverse Compton photons 
is above the threshold for electron-positron pair production 



qE/c ~ {rrieC ) ]. This condition gives 

2 \ i/a 



7c 



2meC' 
3hu)c 



; 3.9 X 10° 



(25) 



where ljc = eB/nieC is the non-rclativistic cyclotron fre- 
quency (e.g., Lovelace 1987). At this Lorentz factor, Ssyn = 
1.3 eV and eic = 2 x 10"eV. 

The acceleration of particles is more complicated in the 
case of an electron-ion jet. The accelerated electrons will 
tend to increase —Jz more than the opposite effect from the 



ions. This will give rise to an inductive ambipolar field. 
It in turn retards the electron acceleration and enhances the 
ion acceleration. However, analysis of this problem is beyond 
the scope of the present work. 



3 CONCLUSIONS 

We have developed a model of Poynting jets as a transmis- 
sion line carrying a net axial current and having a potential 
drop across it. The currents and voltages are of the order of 
3 X 10^* A and 2.7 x 10^" /3 V. The energy of ions accelerated 
across this voltage is large enough to account for ultra high 
energy cosmic rays (Lovelace 1976; Biermann et al. 2001; 
Ostrowski 2002). Further, we derive the transmission line 
impedance and the electromagnetic energy flow in the jet. 
The observed current in 3C 303 is used to independently es- 
timate the electromagnetic energy flow in this magnetically 
dominated jet. Time-dependent but not necessarily small 
perturbations of a Poynting-flux jet - possibly triggered by 
an irregularity in the jet - are described by the "telegra- 
phers' equations," which are wave equations for the current 
and voltage on the line. The voltage and current consist in 
general of forward and backward propagating components. 
The disturbance of a Poynting jet by an irregularity can give 
rise to localized dissipation in the jet which may explain the 
enhanced synchrotron radiation in the knots of the 3C 303 
jet and in the much smaller apparently stationary knot HST- 
1. 

Lastly, we consider relativistic Poynting jets relevant to 
parsec-scale jets such as HST-1. The reflected voltage wave 
from an inductive load (or jet termination) can lead to a 
backward propagating wave which causes the magnetic in- 
sulation to breakdown. That is, it gives |E|/|B| ^ 1. At the 
threshold for breakdown, |E|/|B| = 1, positive and nega- 
tive particles are directly accelerated in the E x B direction 
which is approximately along the jet axis. Particles can be 
accelerated up to Lorentz factors ~ 10^ in a short time in- 
terval of the order of the light travel time across the jet. 
This particle acceleration mechanism is distinct from parti- 
cle acceleration in shock waves and that in magnetic fleld 
reconnection. 

The breakdown of magnetic insulation |E|/|B| ^ 1 is 
not possible in a plasma that is modeled everywhere by 
ideal relativistic magnetohydrodynamics (RMHD). This is 
because the Ohm's law of RMHD E-|-vxB/c = requires 
that |E|/|B| |v|/c < 1 everywhere. Of course an actual 
plasma can readily have |E|/|B| 1. A common example 
is magnetic reconnection where there is a region in which 
the direction of the magnetic fleld reverses so that |B| goes 
through zero on a surface. But even in regions of non-zero |B| 
a plasma can have |E|/jB| ^ 1 transiently as in pulse propa- 
gation along magnetically insulated transmission lines (e.g., 
Shope et al. 1978; Samokhin 2010). The restrictive Ohm's 
law constraint of RMHD can be avoided in relativistic- 
electromagnetic particle-in-cell simulations where the orbits 
of individual particles are calculated (e.g., Lovelace, Gandhi, 
& Romanova 2005). f 

An open question regarding the propagation of relativis- 
tic current-carrying Poynting jets is the kink instability. In 
the context of laboratory current-carrying plasmas, the kink 
instability is predicted and observed to occur if the current 
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is larger than a critical value which is the Kruskal-Shafranov 
condition (e.g., Kadomtsev 1966; Huarte-Espinosa et al. 
2012). The theory of the instability for relativistic Poynt- 
ing jets has not been developed. Observational evidence for 
large current flows in astrophysical jets (e.g., Kronberg et al. 
2011) suggest that the jets can withstand the kink instabil- 
ity. On the other hand a mechanism for the nonlinear stabi- 
lization of the kink instability was proposed by Kadomtsev 
(1966, p. 188). 
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APPENDIX A: WAVEGUIDE MODEL OF 
POYNTING JET 

There are a wide range of a^xisymmetric, translationally- 
symmetric, time-independent Foynting jet models which 
satisfy the relativistic, force-free equation. 



^ ^ 1 d[r\Bl-E',)] _^ 

J.2 fj^j. 



(Al) 



dB 

dr 

(cgs units). This assumes that there is plasma everywhere, 
but that the 'plasma beta' - the ratio of the kinetic energy- 
density to electromagnetic field energy-density is <^ I. The 
time-independence of the solutions of (Al) implies that the 
energy flow in the jet is unidirectional. Here, we assume 
this is the +z direction. Realistic models can in princi- 
ple be derived for example from relativistic particle-in-cell 
simulations (e.g., Lovelace et al. 2005). Here, we discuss a 
waveguide-like jet model alternative to that considered in 
the text. The electric and magnetic fields are similar to the 
fields in a TEM coaxial waveguide. This model is suggested 
by the work of Bisnovatyi-Kogan & Lovelace (1995). 

Equation (Al) is satisfied by the fields shown in Fig. 
Al. For r < ri, the magnetic field is = const and the 
electric field is zero. For n < r < r2, the magnetic field is 
= —Bo{ri/r) < 0, and the electric field is Er = r}B^ < 
where B^ = (1 — rf)^^^'^Bz, and 77 = const < 1 is dimen- 
sionless (cgs units). For r2 < r-z < r:i, the magnetic field 
is B^ = {1 — rf Y^^Bo[r]_/r) and the electric field is zero. 
The condition on the magnetic field at rs is the same as dis- 
cussed in the text. The plasma in the region ri < r < rz has 
a uniform axial E x B drift velocity jycz or Lorentz factor 
F = (1 — rj^)^^^^ . This model involves only a small number 
of parameters, (ri, rg, Bo, '?)■ 

In MKS units, the current carried by the core of the jet 
(in A) is 7o = 27rriBo/Mo- Thus the poloidal flux carried by 
the jet is 

$z = 7rr?B, = i(l 



2n1/2 r 

The voltage drop across the jet is 

fr 

Vo = - 

J n 

where 



(A2) 



drEr 



(A3) 
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Figure Al. Field profiles of a Poynting jet model alternative to 
that of the text (eqn. 2 and Fig. 2). Here, we have taJcen 77 = 0.9, 
ri = 0.5 with a narrow boundary layer, and r2 = 3 also with a 
narrow boundary layer. 



2n' 



-r?ln — 
£0 V^i 



60 r?ln ( ^ ) n [MKS] . (A4) 



This formula is the counterpart of equation (4). The elec- 
tromagnetic power flow in the Poynting jet is 

r-'"2 



_ 27r r 

/^O Jr. 



rdr ErB^ = ZqIq 



which agrees with the equation (3). 

The electric field energy per unit length is 

WE = ^2n f ' rdrE^ = \cVo , 



where 



C = 



27reo 



ln(r2/ri) 

is the capacitance per unit length of the jet. 

The magnetic energy per unit length of the jet is 



(A5) 



(A6) 



(A7) 



Wb 



— / rdrBl + — / 

MO Jo Mo Jri 

Carrying out the integrals we find 

„2 



rdrBj. 



L = 



Mo 
27r 



ri 



+ (1 - ?7 ) In ' 

■ r-2 



(A8) 



(A9) 



which is the inductance per unit length of the jet. 

From equations (A7) and (A9) we find the wave speed 
on the jet is 



Uip = c [ 1 + 



(l-772)in(r3/r2) 



-1/2 



21n(r2/ri) ln(r2/ri) 
The wave impedance of the transmission line is 



Z = Zo{l + 



1 



V 



21n(r2/ri) 



(l-?7")ln(r3/r2) 
ln(r2/ri) 



1/2 



(AlO) 



(All) 



Equations (AlO) and (All) are the analogues of equations 
(12) and (14). 
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